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We analyze dynamics of the FRW models with global rotation in terms of dynamical system 
methods. We reduce dynamics of these models to the FRW models with some fictitious fluid which 
scales like radiation matter. This fluid mimics dynamically effects of global rotation. The significance 
of the global rotation of the Universe for the resolution of the acceleration and horizon problems 
in cosmology is investigated. It is found that dynamics of the Universe can be reduced to the two- 
dimensional Hamiltonian dynamical system. Then the construction of the Hamiltonian allows for full 
classification of evolution paths. On the phase portraits we find the domains of cosmic acceleration 
for the globally rotating universe as well as the trajectories for which the horizon problem is solved. 
We show that the FRW models with global rotation are structurally stable. This proves that the 
universe acceleration is due to the global rotation. It is also shown how global rotation gives a natural 
explanation of the empirical relation between angular momentum for clusters and superclusters of 
galaxies. The relation J ~ M 2 is obtained as a consequence of self similarity invariance of the 
dynamics of the FRW model with global rotation. In derivation of this relation we use the Lie group 
of symmetry analysis of differential equation. 



I. INTRODUCTION 



Rotation is a very universal phenomenon in nature. We can observe rotating objects at all scales of the Universe, 
from the elementary particles to planets, stars and galaxies. The question is, whether this property is an attribute of 
the whole universe at a very large scale structure. On the other hand, if our universe does not rotate, then the question 
is why and how does this happen? Since rotation is generic in the universe, the possible rotation of the universe cannot 
be excluded at the very beginning. Moreover, we should explain the physical mechanism which prevents universal 
rotation of the universe . 

Of course, if the Universe rotate, then there should exist some observational manifestations of this global rotation. 
Following the standard approach we can study the motion of test particles and photons on the background of cor- 
responding spacetimes which admit the rotation. Our idea is to test the effects of rotation by classical cosmological 
tests like luminosity distance or angular size of radiogalaxies. 

Since the work of Lanczos , Gamov Q , Goedel Q , and Hawking || , the cosmological models with rotation have 
been studied as well as the behavior of geodesies in such spacetimes. Although the quite strong upper limits for 
cosmic vorticity were obtained from analysis of CMB or BBN 0, 0, 0> @ ai l these works based on the model in 
which shear and vorticity effects are inseparable, i.e., in the sense that zero shear automatically implies zero vorticity. 
Therefore, the obtained limits actually placed not on the vorticity, but rather on the shear induced by it within the 
specific geometry of spacetime . Finally, one thus needs a separate analysis of these limits in expanding cosmological 
models with trivial shear (shear- free) but non-zero rotation and expansion (it is an idea of Obukhov pj). Principally, 
there are two classes of such models in which this analysis can be performed 

1. Newtonian shear- free models so-called Heckmann-Schiicking models [Tol fTT| : 

2. general relativistic spatially homogeneous models with the geometry of the Bianchi IX and the cosmological 
constant 0. 

Let us consider the first class of models. We have a homogeneous Newtonian universe defined on product of the 
three-dimensional Euclidean space and absolute time coordinate. This universe is homogeneous, density and pres sure 
of the fluid have no spatial dependence and the velocity vector field depends linearly on spatial coordinate |l2l Il3| . 
In this case contrary to the general relativity theory there exist shear-free solutions which satisfy the Heckmann- 
Schiicking equations and describe both expanding and rotating universes [l^. The rotating homogeneous non-tilted 
universes when considered on a relativistic level filled with perfect fluid must have non- vanishing shear |l5l Il6| . 
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Therefore, it is very difficult to consider both effects of anisotropy and rotation, that's why we study the Newtonian 
counterpart of model with rotation. Moreover, it is reasonable to assume that the shear scalar is sufficiently small 
compared with angular velocity scalar since the shear falls off more rapidly than the rotation as the universe expands 

To compare the results of analysis of supernovae data with a general relativistic model we formally consider a 
non-flat model Qk,o 0, although the satisfactory interpretation of £lk,o can be found beyond Newtonian cosmology 
in the general relativity. 

The motion of the fluid in a homogeneous Newtonian universe is described by the scalar expansion 8, the rotation 
tensor u) a b, & n d the shear tensor a a \,. The homogeneous rotation of fluid as a whole is usually called the global rotation 
of the universe |Ts| . 

The problem of the rotation of the whole Universe attracted the attention of several scientists. It was shown, that 
the reported values are too big, when compared with the CMB anisotropy. Silk [24| pointed out that presently the 
dynamical effects of a general rotation of the Universe are unimportant, contrary to the the early universe, when 
angular velocity loq > ICU 13 rad/yr. He stressed that now the period of rotation must be greater than the Hubble 
time, which is consequence of the CMB isotropy. Barrow Q also addressed this question and showed that the cosmic 
vorticity depends strongly on the cosmological models and assumptions connected with linearisation of homogeneous, 
anisotropic cosmological models over the isotropic Friedmann universe. For the flat universe they obtained the value 
(u/Hq) = 2 • I(T 5 . 

While the investigations of the cosmological models with global rotation has a long history, the global rotation 
of the Universe is still not detectable by observation. Some observational evidences for the global rotation of the 
Universe have not been confirmed [20l l2lL 12*2^. The cosmological origin of the rotation of astronomical objects was 
analysed by Fil'chenkov in the framework of quantum cosmology and tunnelling approach |23| . 

The alignment of galaxies along the direction of the global rotation would also be some indication of the discussed 
phenomenon. However, the primordial dipole anisotropy of the distribution of rotation axes of galaxies is blurred due 
to the irregular shape of the protogalaxy, which can lead, after phase transition, to an almost random distribution. It 
is clear that the angular momentum of a structure is a sum of spins of components and its own angular momentum. 
A lot of work has been done in study of the alignment of galaxies in superclusters (for review see |25j). Here, the 
result are ambiguous, but several independent investigations claimed the existence of the weak alignment of galaxies 
in respect to the supercluster main plane [26l l27l I2H , |29| . 

I ndep endent observations of Supernovae type la (SNIa) have indicated that our Universe presently accelerates 

Usui mm. There 

is a fundamental problem for theoretical physics to explain the origin of this acceleration. One 
of the first explanations was the cosmological constant. It is also possible to introduce the global rotation as we do 
in the present paper. While such a possibility seems to be attractive in the context of the accelerating universe, the 
cosmological constant is still needed to explain SNIa data because the contribution to acceleration coming from the 
global rotation is too small |34| . 

In our previous papers we study some constraints on rotation coming from SNIa data, primordial nucleosynthesis as 
well as from the age of the oldest in the Universe |34j . The main aims of this paper is rather theoretical than empirical. 
We show that the analysis of the model dynamics is more effective after rewriting the basic equation to the new form 
with dimensionless quantities. Then we obtain a two-dimensional dynamical system in which coefficients have a 
simple interpretation of a density parameter of fluid fulfilling the Universe in the present epoch. In this approach 
dynamics is described by the Friedman-Robertson- Walker (FRW) model with some additional noninteracting fluid 
for which we can define density parameter fl^ = p^jZH^. The presented formalism gives us also a natural base 
to discuss the influence of global rotation on acceleration of the Universe. From the observational constrains it is 
possible to find a limit for the angular velocity of the Universe. For example, Ciufolini and Wheeler [35| obtained 
ujq — 6 • 10~ 21 rad s _1 ~ 2 ■ 10~ 13 rad yr from the CMB limit for the homogeneous B(IX) models with rotation. The 
limits for angular velocity of the Universe obtained from the CMB and SNIa were also discussed j3^|. Here we also 
investigate some dynamical effects of the global rotation on the solving classical problems in cosmology, namely the 
horizon, flatness and cosmological constant problems. 

Organization of our paper is as follows. In section 2 we present Hamiltonian dynamics of the FRW model with 
global rotation. In section 3 we discuss some properties of exact solutions with rotation. In section 4 we discuss the 
domains of the phase plane for which the cosmic acceleration and horizon problem are solved. In section 5 we construct 
the phase portraits for the FRW model with rotation and dust (7 = 1). We also investigate general properties of 
dynamics by showing that generally it can be reduced to the form of the FRW models with non-interacting multifluid 
scaling like radiation. In section 6 a detailed analysis of dynamics on phase portrait is given in two dimensional 
phase plane. In section 7 we discuss the flatness and cosmological constant problems in the context of model with 
rotation. In section 8, it is proved (without of any help of empirical evidence) that formulae J ~ M 2 is exactly 
valid for the universe with dust (7 = 1). For proof we use the Lie group theory of symmetry analysis of differential 
equations, describing the evolution of the universe. Then the relation between angular momentum and mass becomes 
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a constraint relation between invariants of a self similarity group of symmetry of the basic equations. In section 9 we 
summarise obtained results. 



II. REDUCED HAMILTONIAN DYNAMICS OF FRW UNIVERSE WITH GLOBAL ROTATION 

In this paper we adopt the Hamiltonian formalism. It gives immediate insight into possible motion because the 
dynamical problem is analogous to that of a particle of unit mass moving in the one-dimension potential. 

The Raychaudhuri equation describing the relation among scalar expansion 0, the rotation tensor uj a b and shear 
tensor a a h_(jf the acceleration vector A a = u b W v u a vanishes, which would necessarily follow in the case of dust) has 
the form |35l | 

6 = -\<d 2 - 2(ct 2 - oj 2 ) + A - 47rG(p + 3p) (1) 

where u> 2 — uj a bdJ ab /2, a 2 — a a bC at '/2, A is the cosmological constant, G is the gravitational constant, the fluid is 
perfect fluid with the stress-energy tensor T a b = (p + p)u a Ub + pg a b (p is the mass density and p is its pressure). The 
shear scalar a 2 vanishes for the FRW model, while for the anisotropic Bianchi I or Bianchi V model (hereafter B(I) 
and B(V)) it satisfies the following equation 

^ = -29a 2 . (2) 
dt y ' 

Our aim is to consider the FRW model and its generalisation — anisotropic B(I) or B(V) models with global rotation. 
In the anisotropic generalization we have 

36 = [E x + H 2 + H 3 ), 

where Hi = <i(ln ai)/dt are the Hubble functions in three different main directions, and a, (i = 1, 2, 3) are corresponding 
scale factors. Equation Q can be also obtained within the Newtonian theory. However, in general relativity for the 
perfect fluid model with vanishing shear a — 0, and acceleration u a = 0, the field equations can be consistent only 
if w0 = 0. For the dust matter fluid with a — 0, u a — 0, particles of dust are moving along geodesies. This is in 
contrast to the Newtonian case [l7j . Therefore, expanding and rotating spatially homogeneous universes filled with 
perfect fluid must have a non-interacting shear. Li |l8| considered the case that a is sufficiently small compared with 
lo since the shear following equation © falls off more rapidly than the rotation, during the expansion of the Universe 
0. In the generic case in general relativity there is no solution with rotation and without shear. In other words, it has 
been shown that spatial homogeneous, rotating, and expanding universes with the perfect fluid have the nonvanishing 
shear [l7ll35l|. However, there exists such a possibility in the Newtonian cosmology because some standard obstacles 
of Newtonian cosmology can be overcome in this model . The conservation of energy gives in a relativistic case 

P = -e(p+p). (3) 

The occurrence of the term p in the factor (p + p) has a special relativistic effect. It is a consequence of the inertia 
assigned to all forms of energy. The Newtonian counterpart of equation Q corresponds to p — (dust for which 
world lines of fluid are geodesies 01)- 

For a rotating, spherically symmetric, self-gravitating system with a characteristic size R, mean density p, and 
rotation uj, the angular momentum J is given by J = 2MR 2 u>/5. Therefore, the angular momentum conservation 
can be written in the Newtonian case in the form 

uipa 5 = const (4) 

where a 3 = aia2&3 is an average scale factor. Principle (0} has its counterpart in the relativistic case 

uj(p + p)a 5 = const. (5) 

In the general relativistic case rotation is governed by a pair of nonlinear propagation and constraint equations which 
for the case of the barotropic equation of state reduces to |n| 

!(i-|die-. 
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where uo a is a shear eigenvector (<J a b^ 2 = cu; a ), and c s is a velocity of sound (c 2 = 0). 
The Einstein equation for the considered case has a first integral in the form 

- 2 6 2 8ttG 2 2 ct 2 3 R A 

* = t = — ^-r +T-1T+3 (6) 

where 3 i? is the Ricci scalar of three space equalled to 6k/a 2 (k = —1 for B(V) and k = for B(I)) and for a natural 
system of units 8nG = c = 1 . We usually work with the equation of state of the form 

P = (7 - 1)P (7) 

where 7 is a constant. From the physical point of view, the most important values are 7 = 1 (dust), 7 = 4/3 
(radiation), 7 = 2/3 (strings), 7a = 27 (brane effects), 7 = 1/3 (solid dark energy or topological defects). Let us note 
that in equation © there is no dependence on a particular form of the equation state. 

If the matter is in a pressure-free form it is represented only by its 4- velocity u a and its energy density p. From the 
momentum conservation we obtain that matter moves along geodesies. For our further analysis it would be useful to 
rewrite equation to the form of a differential equation 

6j = -2ujH or u = (37 - 5)ujH (8) 

for the equation of state in the form J7J. After differentiation both sides of equation (jBJ with respect to the time 
variable and then the substitution of equation © we can obtain the Raychaudhuri equation. It explicitly demonstrates 
that in terms of the average Hubble function 

H =^(H 1 +H 2 + H 3 ), 

equation (JSJ) is really a first integral of basic equation Q . In the considered case we write this equation in the form 

F = -H 2 -^a 2 - W 2 ) + |-ip(3 7 -2), (9) 

where the bar over H will be further omitted. For the FRW symmetry we have Hi = H2 = H3, (but in general the 
Hubble function H is understood to be averaged over three different main directions). For anisotropic models B(I) 
and B(V) the average Hubble parameter H is calculated from the average scale factor H = \n(d/dt( &a\ «2 ^3) ) • Hence 
equation 10 for the model with the FRW symmetry can be rewritten to the form 

- = -V-- 2 ) + ^p(37-2) (10) 
a A 3 6 



where 



Poa 37 , a = oo a 3 > w = ujqci 5+37 . 



Let us note that from equation J8J) we obtain that u>(a) does not preserve its form in different epochs. The value of 
J depends on 7 ( J = topja 5 ). Of course 

i : = u L + p + 5H = ±_ 3iH + 5H = 0: 

J LU p LU 

and conservation of J takes place if only lu obeys relation JSJ. 

Therefore the right-hand side of (|10fl can be expressed in terms of the scale factor a (or average for B(I) and B(V)) 

.jjofc-e _ 5l_L^ poa -37 + A + ^-ro+e^ a (n) 

Equation 111|) can be rewritten in the form analogous to the Newton equation of motion in the one-dimensional 
configuration space {a: a 6 R+} where a is the scale factor 

9V 
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where the potential function for 7 7^ 4/3 is in the form 



v W - -a„-« _ m a -3-, +2 _ _ I^l^o + ^ (13 , 



while for 7 = 4/3 it takes the special form 



7 -4 PO -2 A 2 



7(a) = --^a-* - ^a-' - -a* -Clna + V (14) 

6 

where Vq = const and C is an undetermined integration constant in the Hamiltonian which the energy level should 
be consistent with the form of the first integral ©. Equation (|12|l has a first integral which can be expressed for 
7 ^ 4/3 in terms of the potential function V(a) as 

V(a) + £=V a -± (15) 

Now we can construct the Hamiltonian 

H=^ + V(a) (16) 

Trajectories of system 111 61) are situated on the distinguished energy level, Ti = E = const. We can also define a new 
Hamiltonian flow in the one-dimensional configuration space on the zero energy level by translating the curvature 
contribution to the Hamiltonian. Then 



n=- + V(a)=Q (17) 



where for 7 ^ 4/3 



and for 7 = 4/3 the Hamiltonian takes the special form 

d 2 , , (7 2 a On , A , „ k a 1 

W = T +V(«) = -f« -6° >-Cl*a + - + Y = E=- r > 

which becomes consistent with the form of first integral JSJ only for C = 0. Therefore, for 7 / 4/3 the physical 
trajectories lie on the zero-energy level Ti = E = which coincide with the form of the first integral Q18[). For the 
special case of 7 = 4/3 trajectories lie on the distinguished energy level Ti — — w 2 /3 = const < 0. In both cases 
dynamics is given in the Hamiltonian form and trajectories lie on the distinguished level of constant energy. Finally, 
for radiation (7 = 4/3) we obtain the new Hamiltonian considered on the zero energy level 

oo -4 Po -2 , a 2 A 2 (k 1 2 \ , . 

Ti = -a a H a + — — lot, = 0. (19 

6 6 2 6 V 2 3 V V; 

Note that this model is dynamically equivalent to standard model with some effective curvature, 

k cS = k + ^iv 2 . (20) 

Let us consider the FRW (o~n = 0) model filled with some pure unknown matter X such that px = Wx Px ■ Then we 
obtain 

V(a) = -f a -^ 2 + ±a 2 - ^ a 2 ^^ + * (21) 
6 6 iyswx — 4) 2 

The mixture of noninteracting dust and unknown matter X satisfies the equation of state 

P = Px + = (7(a) - 1)P = ggo ffl3( r*-i) T I Pi P = Pmoa- 3 + P xoa- 3wx . (22) 

Pxo 
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Then the potential function is given by the integral 
V(a) = 



3(7(0) -2) 2 2 , . A 2 2 _ 10 / fj(a) , 
- AZi -£ " + 3^ (a) - g - 3^o« exp 6 / -^da 



a da 



(23) 



where 7(a) is derived from equation 12 I'll and for the mixture of dust and unknown matter we can find that 



[M da= 

J a J eoxa 



-da 



w x 
m + 3 



ln(e 0m + e ox a m+i ) 



(24) 



where m = —3(1 + wx)- Finally we obtain 



V{a) = _PXl a -^ X+ 2 _ ^ fl -4 _ A 2 k 2 

W 6 6 6 2 3 °J (e 0m + eoxa m + 3 ) 2 



(25) 



Let us consider for example the case of solid dark energy (wx — 1/3, m = —4). To do this we substitute these 
values into integral (|24|l which gives 



J = 



da 



1 



1 



{-Ifn 



fc„,fc-i 



6 , 1 + na 



a 8 (e 0m + eax) 2 eom£ox a 8 (l + cra) (7 - k)\a 



'In- 



where a = e m/eox- 

Then we obtain potential in term of the integral J 



V^(a) = -- — a -a 4 a H wSJ. 

w 6 6 6 2 3 



(26) 



III. EXACT SOLUTION WITH ROTATION 

In the case of the dust matter, an exact solution with rotation seems to be especially interesting. It can be easily 
shown that it can be simply obtained if we note that the corresponding dynamical equation reduces exactly to the 
best known case of the FRW model with dust, radiation, and the cosmological constant. 

In the case of vanishing shear and the cosmological constant A = an exact solution can be given in the simplest 
form. There are three possibilities of positive (k — +1), zero (k = 0), and negative (k = —1) curvature. 

a) If k = +1 equation l|17l) becomes 



-2V(a) 



Pa 1 2 loI 
3 a 3 -r~ 2 3 (37 - 4)a" 6 T+8 



-1, 



and since d 2 has to be positive there are maximum and minimum values for a (if 7 < 4/3) as can be derived using 
elementary techniques. The general solution for 7=1 (dust) is 



,.,/, a . a — 2a(t) _ 
t = -(j3 + aa- a 2 ) 1 ' 2 - - arcsm - — | - j ' + C, 



4/3)V2 



where n = p /3, (3 = —%Wq and a min < a < a max = \{(n 2 + 4/3) 1/2 + a), lu < p a V2. 

Since only the first derivative vanishes for a = a max , we obtain an oscillating singularity-free universe, 
b) If k = equation l|17|) gives 

■ 2 a P 



7,37-2 



-67+8 ' 



Here d 2 is always positive for a > —/3/n. So there is no maximum value for a. The result is a monotonically expanding 
universe without an initial singularity. 
For 7=1 the general solution is 



t 



:(na + l3) 1/2 (na-2p) + C. 
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c) The last case, k = — 1, gives the equation 

•2 = a P _i_ 1 

° a 3 "/- 2 a- 6 ^+ 8 ■ 

This results is also a monotonically expanding universe without an initial singularity. 
The scale factor has a minimum 



_ -a + ^a 2 - 4(3 



For 7 = 1 the solution of i|17|) is 



if a 2 - 4(3 > 0, 



if a 2 -4(3 = 0, 



t= (a 2 +aa + (3) 1/2 - ^ In (2(a 2 + aa + (3) + 2a + a) + C, 



t=(a 2 + cta + (3) 1/2 - - In (2a + a) + C, 



t=(a 2 + cta + /?) 1/2 - - In (2a + a) + C. 



Exact formulas for observables are also possible. For example, in the case of A = the luminosity distance Dq 
between a galaxy, which emitted light at the cosmic moment t and the observer who detected the light at the present 
moment t , is given by 

_ n m , z + (n m , + 2q,,o - 1) {[2iy z 2 + 2z (2q,, + n m , ) + 1] 1/2 - 1} 

^o[^m,o + 2^w,o (2^,0 + 2Sl m> o — 1)] 

The above formula is valid for all values of curvature index (see for details 13). 

In the more general case of the nonvanishing cosmological constant we can find also a solution in the term of an 
elliptic function. In the case of dust matter, the problem of finding an exact solution is equivalent to integration of 
the equation 

^ -y 2 2 



which takes the following form in the conformal time 



= ^u,0 + £l m fiX + ilkfiX 2 + flAfiX 4 = W{x), 

where dr/ = dt/x, x = a/a and ao is a present value of the scalar factor a. 
Then it can be integrated explicitly using the Weierstrass p function, to yield 

_ y/W{x)W{x ) + W(x ) 4lxl + 2kx a + m 6lxl + k 
P( - V) ~ 2(a;-xo) 2 + 4{x - x ) + 12 

or equivalently 

y/Wfajp'(r)) + i[p(ri) - 6 " A '° j> +n "- ] (4O A , ^ + 2Q M a;o + n ro , ) + W{x )n A , x 



X — Xq 



2 [ p ( v ) - 6"A,o» +" fc .° ]2 _ ln A . W(x ) 



where xo = x(i] = 0), and p is constructed with the invariants gi = ^A,o^u,o + j^^fcO' 93 = |^A,o^fc,o^w,o — 
(fi fe ,o/6) 3 - ^OA,0^m,0- 

In the considered case the effect of the global rotation is modelled by negative radiation. Then, if we put in the place 
of dimensionless fl r .o some effective constant r^ ff , (which combine all contribution coming from matter scaling like 
radiation) we immediately obtain the exact solution in terms of the Weierstrass elliptic function [3^ . It is interesting 
that in this case formulas determining the magnitude-redshift relation can be also given in the exact form. After 
introducing k c ff (see formula Q2U[ll we can also simply obtain the corresponding exact formulas from that obtained for 
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IV. COSMIC ACCELERATION AND HORIZON IN THE UNIVERSE WITH GLOBAL ROTATION 



In our further analysis we will explore dynamics given by the canonical equations. Assuming x = a and y = a = x 
we have from equation l|17|) 



y 



&H _ 
dy 

m 

dx 



!J 



dV 
dx 



(27a) 
(27b) 



with the first integral in the form y 2 /2 + V(x) = 0. Then, we can perform a qualitative analysis of autonomous system 
(I27|l in the phase plane (a, d) = (x, y). 

The general observation is that if we consider the two-dimensional Hamiltonian dynamical system then eigenvalues 



of the linearization matrix Ai, A 2 : X — — ^r\ ao - They are saddle points if ^r|a < and then eigenvalues are real of 
opposite signs otherwise they are centres at which eigenvalues are purely imaginary and conjugated. This dynamical 
system rewritten to more useful form will be considered in section 5. 

We can also observe that trajectories of system Q27JI can be integrable in quadratures, namely, from the Hamiltonian 
constraint H. = E = 0, we obtain 



t-to 



da 



(28) 



Note that it is possible to make the classification of qualitative evolution paths by analyzing the characteristic curve 
which represents the boundary of domain admissible for motion. For this purpose we consider the equation of zero 
velocity, a — which represents the boundary in the configuration space. Because 



2V(a) 



(29) 



the motion of the system is limited to the region {a: V(a) < 0}. Let us consider the boundary set of the admissible 
for motion configuration space given by a condition 



dM = {a £ R+ : V(a) = 0}. 



(30) 



For the special case of 7 = 4/3 of course levels V(a) = E should be considered. From the exact form of the potential 
function (|18|l or (|23[l parameters A or ujo can be expressed as a function of a. For example from (|18|l for 7 ^ 4/3 we 
obtain 



A(a) = ^ 



6 6 



1 ^0 -2(37-4) , £ 

3(37-4) 2. 



and for 7 = 4/3 



A(o) = 



7 -4 PO _2 



-~6 a 



6 



(31) 



(32) 



Finally, we consider the evolution path as a level of A = const and then we classify all evolution scenarios modulo with 
their quantitative properties of dynamics (to compare see |39L l40j | ) . For the special case of radiation the corresponding 
classification of the FRW model effective curvature k e g can be achieved in a similar way. Then we obtain the analogous 
phase portraits is shown on Fig. 

From the physical point of view it is interesting to answer the questions: are the trajectories distributed in the phase 
space in such a way that critical points are typical or exceptional? How are trajectories with interesting properties 
distributed? For example, along which trajectories the acceleration condition, a — —dV/da > is satisfied? In order 
to fulfil that V(o) must be decreasing. One can easily observe this phenomenon from the geometry of the potential 
function. In the phase space, the area of accelerations is determined by condition y > or by corresponding condition 
in the configuration space 



3 7 -2 
6" 



-p a 



-37 



A 
3" 



-10+67 



< 0. 



(33) 



The presence of the last term in inequality l|33|) demonstrates that in the configuration space the domain of acceleration 
is larger than for the case of the FRW model with vanishing global rotation. Let us note that the positive value of 
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the cosmological constant A acts in the same direction as rotation, i.e., global rotation acts as dark energy (or dark 
radiation). 

It can easily be demonstrated for the FRW models that if a(t) — > const as a — > then the model has no particle 
horizons (in the past) |4l| . Indeed, if there exists a constant C such as for sufficiently large p, the velocity of the scale 
factor is upper bounded da/dt <C and 

a to 

[*<C[£. (34) 
J a J a 

o o 

The integral on the left hand-side of Ij34(l diverges and there are no causally disconnected regions. Putting this in terms 
of (x, y) variables one needs the condition, as x — * 0, y — ► const, as the sufficient condition for solving the horizon 
problem. Now, from the Hamiltonian constraint we obtain that as x — > 0, then V(x) — ► const (zero is included). 



V. THE PHASE PLANE ANALYSIS— GENERAL PROPERTIES OF THE MODELS 



Let us apply a dynamical system method to analyse the system under consideration. The dynamics of cosmological 
models with global rotation, the cosmological constant A and matter, satisfying the equation of state 0, is described 
by the following dynamical system 

a = y (35a) 

y = - 2 r la-* - h_l poa -^ + A fl + 2^6,-9. (35b) 

Let us note that some general properties of system <|35[) for vanishing shear (<7q = 0) and the matter in the form of 
solid dark energy (or topological defects) with 7 = 1/3. Then when a — *■ the effect of the cosmological constant is 
negligible. In this case y = ^po + |wga -7 and the term with rotation cannot dominate the matter term with positive 
energy because V(a) < 0. From the existence of the first integral we have 



y = y/-2V(a) = J ^p a + ^wgo" 



and a > ( - — uj 



2 ^ 1/7 



2 

3po 



It is possible that density of solid dark energy and the rotation compensate each other, when the scale factor reaches 
a = flmin near the singularity. 

As a — > we have y — ^po, i.e., y ~ t and a ~ t 2 and we obtain that the horizon problem is solved because 
y = a — » as t — > 0. The phase portrait for this case show Fig. 

To apply the method of dynamical system let us consider the FRW model with dust matter, the cosmological 
constant, and global rotation. Then the corresponding dynamical system has the form 

a = y (36a) 

y = -ip a~ 2 + g-a + ^a -3 (36b) 

with the first integral, for the general form of the equation of state p = (7 — l)p, given by 

2 1 -1 A 2 2 2 a 2 ^- 4 > 

V =~ 3 Poa +-a (37) 

Near the initial singularity of type, i.e., as a — ► const the generic flat solution without the cosmological constant 
has the form 



If we assume that ujq <C 1 then from solution (|38|l we obtain a simple monotonic solution with the singularity at 



an 

Pa 




t-to = — — Wa + — Si. (39) 



10 



Note, that as in the previous case there is the possibility that the scale factor reaches a — a m i n in the initial 
singularity. From the condition that motion is admissible in the region V(x) < we have for the flat model the 
condition 

1 i 1 _, 

\fx: - -tt m , x - -ItvflX <0 

or in terms of redshift z 

Vz: n ro ,o + fi W)0 (l + z) > °- 
For the minimum a = a min and the maximum p m we finally obtain the bound 



1 + Z. 



SJllg 



where O mi o — 0.3 and z s ing is the value of redshift at the moment of singularity a = a m ; n such that 1 + z s - ms = 

It would be useful to compare the present effects of radiation and rotation because both scale in the same way. The 

rotation effect could dynamically dominate radiation i.e., 

-n u , > n rfi ~ 0.5 ■ io~ 4 . 

The contribution coming from both radiation and rotation we denote as f2*°Q al . Then, from equation l|36b(l in the 
case without the cosmological constant it can be obtained the relation 

x = -— x [ 1 + ^7 X )■ 

Therefore, for xq — — n 7,, ° we have an inflection point on diagram x(t) — which separates domains in which x 
has different convexity. 

Our starting point for further analysis of results will be dynamical equations rewritten to the new form using 
dimensionless quantities: 

x=-, T=\H \t, fi, i0 =*!°-, (40) 

AO Pcr,0 

P = Po — , o- = o-q — , U) = Uq[ — ) , (41) 

\ «o / \ao/ \ao/ 

with H = p cr = 3Hq 2 where the index denotes present day values (at time to). After introducing the density 
parameter of rotation for 7 7^ 4/3 the dynamical effect of global rotation is equivalent to some kind of additional, 
noninteracting fluid for which we have 

Puj = ^L^(37-5) = ^^-10^ ^ = P»_ = n ^ a .6 7 -10 j (42) 

37-4 p cr 

10-67 , n 

= , p u = {w u - l)p u , (43) 

and negative energy and in consequence the negative density parameter f2 w . Hence, the basic dynamical equations 
are 



l^,o + \ E n ifi x 2 - 3w > (44) 



2 

»=i 



i; = x (45a) 
y = i = ^^,o(2-3«; l ) a ; 1 - 3 ^ (45b) 



11 



where for Wi 7^ 2/3 (7 7^ 4/3) 

v(x) = -in*, - ~ ^ "i.o^" 3 ' 01 . ( 46 ) 



2 

i=l 

lo - 2 go 

2 3(3 7 -4) : 

and the system is now defined on the zero energy level Ti = E = 0. 
If = 2/3 (7 = 4/3) we have 



V{x) = \ ^o^-^ + ^o 2 , (47) 

where in general px = (wx ~ & n d for unknown dark energy X, w\ = 1/3 (solid dark energy), W2 — 2/3 (effects 
of string), u> 3 = 7 (perfect fluid), wa = (cosmological constant), = 3 7 (for rotation), w>, = 27 (effects of 
brane). 

We could in principle include the curvature term ^Qk,o into the sum in the right hand sides of <|44f) . it would 
correspond to Wk = 2/3. We see that in particular the case xq = 1, xq = 1 is a solution. From this definition 
xq = — go, where qo is the deceleration parameter. Therefore, J^^i.o + ^fe,o = 1 (presently but this equality can be 
extended to any given time). For the special case of the universe filled with radiation matter we can also formally define 
energy density of some fictitious fluid which mimics the rotation effect and scales like "curvature fluid" , p u = — 2uj^a~ 2 . 
Hence, the corresponding density parameter is 

ft- = P " 



and because the system is defined on the levels 7i — E — we have the constraint i\o = 1. 

Therefore in any case the system under consideration can be represented in the form of the multifluid FRW dynamics. 
Of course, dynamical system (|45|l is the Hamiltonian 

n= P -f + V{x) (48) 

where p x is the momentum conjugated with the generalised coordinate x. 

In the following we consider a universe filled with dust and some unknown component X with negative pressure 
called "dark energy" , so that the strong energy condition is violated px + %Px < 0. We assume that the weak energy 
condition px + px > is satisfied, so wx > for the hypothetical quintessence matter. Hence, the Universe is 
accelerating provided that the following condition is satisfied 

2a(x) = f2 x ,o(2 - Sw x )x 1 - 3wx - Q^qx' 2 + 2Q ufi (3S(x) - 4)x 3(27 ~ 3) + 2fl Afl x > 0, (49) 

where 

S( x ) = n n 7, 7' P = S(x)p 



n 



X,0 



and p, p are total pressure and energy density; p = + wxpx, P = pm + Px- 
Our universe accelerates at present (x = 1) if 

Ox,o(2 - 3wx) - a„,o + 2^,0(3(5(1) - 4) + 2tt A ,o > 0. (50) 

Let put Qa,o — and matter X denotes a candidate for the dark energy (the pure cosmological constant corresponds 
to wx = — 1). 

Let a = ^m,o/^x,cb 8 = ^u,o/^x,o, then is required 

(86 + a - 2)(a + 1) 
Wx < 3(26-, -1) ' (51) 



If a = 3/7 then for small 6 we have 



wx < ^(l + 40.6|<y|), 6<0. (52) 
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Therefore, the corresponding conditions of negativeness of coefficient of state for dark energy is weaker then for the 
case of vanishing rotation. Present experimental estimates are based on estimation of barions in clusters of galaxies 
giving Q m fi ~ 0.3 |42|. On the other hand, the location of the first acoustic peak in the CMB detected by the 
experiments Boomerang, Maxima and WMAP suggest a nearly flat universe. It also implies that our Universe is 
presently accelerating for a wide range of values of coefficients of the equation of state, roughly wx < 0.5. Let us note 
that now we have the constraint on the angular velocity 



;3 (fife o + ^m,o + ^A,o — 1) TT —3Q u ,o TT , r n\ 
= y ^ H ° = V ~ ( - 53 ' ) 

which is obtained as a consequence of the constraint relation = 1- 

It is easy to obtain that for fl Ut o ~ 0.01 which corresponds to the limit obtained by us from SNIa data and the 
present value po = 0.8 • 10 -29 g cm" 3 , luq ~ 2.6 • 10~ 19 rad s _1 that is in a good agreement with the observational 
limit HfJ. 

Let us now make some important remarks concerning acceleration regions transition in the configuration space 
admissible for motion x m i n < X < oo. The "physical region" is x m in < x < 1 because x = a/ao. Accelerated 
expansion starts from same x a given by non-zero solutions of the equation 

dV(x) 

h(x) = ^T 1 ^ = -tt m ,ax + 2Q A ,ox 4 - 20^,0 = 

which represents the boundary of acceleration region h(x) > in the configuration space. 

The function h(x) satisfies the condition h"(x) > if only SIa.o > and for large and small x where it has 
asymptotic forms h(x) oc £l\,oX 4 and h(x) = — £l m fiX — 2^1^.0, respectively. Additionally, we have h(0) = — 2f2 w .o > 0. 
Therefore, h(x) has a minimum at x min = y/^l m ,o/^^A,o which coincides with the same value for the standard case 
with vanishing f^.o- 

Let us consider the model without rotation f2 Wj o = 0. Then x G [0, 1]. Now, we can simply define the fraction of 
acceleration during the whole evolution up today as 



1 - Xa , fl 

p = 1 = 1 _ « ■' 



where x a = ^ is the value of the relative scale factor in the transition epoch a tr . For these models accelerated 
expansion starts at x a , which corresponds to transition redshift 



z a — x a X 1 — \ ~P> ~ ^ 

The fact that z a is so close to zero (for the FRW model with the cosmological constant term, z a — 0.414) is the cosmic 
coincidence problem — z a is shifted towards smaller redshift. 

Let us note that this obstacle is less dramatic in the class of models with rotation f2 Wj o ^ 0, because due to presence 
of rotation contribution z a is shifted towards larger redshift. 



VI. DISCUSSION AND DETAILED ANALYSIS OF THE PHASE PLANE 



Dynamical system methods of analysis dynamics of the FRW model with global rotation seems to be useful if we are 
interested in the qualitative properties of dynamics and influence of global rotation on the dynamics of the universe. 
In this approach we deal with the full global dynamics of the universe whose asymptotic states are represented by 
critical points of the systems. In such a representation the phase diagrams in a two-dimensional phase space allow 
us to analyse the acceleration in a clear and natural way. It is the consequence of representation of dynamics as a 
one-dimensional Hamiltonian flow. 

From the theoretical point view it is important to know how large the class of accelerated models is, or what kind of 
qualitative behaviour of trajectories introduce the global rotation. We will call this class of accelerated models typical 
(or generic), if the domain of acceleration in the phase space driven by effect of global rotation is non-zero measure. 
On the other hand, if only nongeneric (or zero measure) trajectories are represented by accelerated universes, then 
the mechanism which drives these trajectories to accelerate should be called ineffective. Such a point of view is a 
consequence of the fact that, if the acceleration is an attribute of a trajectory which starts with a given initial condition 
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it should also be an attribute of trajectories which start from nearby initial conditions. From recent astronomical 
measurements we obtain that while global rotation effects give rise to the acceleration of the universe, the cosmological 
constant is still needed to explain its rate. 
The system under consideration is 



x = y (54a) 

dV 
dx 



V = ~\ X>,o(3«>< - 2)x 1 - 3 ^ =~ (54b) 



with first integral 



In general, it represents the dynamics of FRW models filled with noninteracting multifluid for each i component in 
the equation of state pi — (w; — l)pi is satisfied. 

For the FRW model with global rotation, the cosmological constant and dust we have 

IO-67 



and 



Woj = jj , w k = -, 7 



x = y (55a) 
y = -^m,oX~ 2 - &A,oX - Qu>,ox~ 3 (55b) 
System (|55|) has the first integral given by 

y 2 = --fim.oz -1 + ^fe,o + Qa,qx 2 + il^^x^ 2 . (56) 

The classification of all admissible evolutional paths in the configuration space is demonstrated in Fig. The phase 
portrait which visualizes the all evolutions for all initial conditions is presented on Fig. [21 For comparison the phase 
portrait with the cosmological constant, and topological defects is in Fig.^ The concordance ACDM model is shown 
on Fig. H 

The critical points of l|55|) are solutions of the equations 

v(x , y ) = o n^.ox- 2 + n^ox- 1 + n k . + n A ,ox 2 = o (57) 

and 



dV 

dx 



= -Q u , x 3 - -Q, m , x 2 + n Afi x = (58) 

(xo,3/o) 



where ^ Oj,o = 1. 

Trajectories of the system belong to the domain admissible for motion 

V = {{x, y) e R 2 : x min < x and V{x) < 0}. 

Therefore, critical points lie on the boundary dT> = {x: V(x) = 0} 

V(x) = ~\ x ~ 2 (Clu,0 + flmfiX + Qk$X 2 + Qa^X 4 ). 

To find the solution xo of (|57|l and l|58|) we consider the probe equation which helps us to formulate necessary 
condition for xq ^ to be solution if only 

-fi w ,o + flkfiX 2 + 3fiA,o^ 4 = 0. 
Of course, if both 17a, o an< i ^fc,o ar e positive then there is no critical point in the finite domain. 
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Let us consider the case of Qk o < 0. Then 



x 



\ 



?yo = o. 



60 A ,o 

where it is assumed that rotation is sufficiently small to ensure that real xq exists, i.e., 



n 2 

120a n 



In the opposite case if the contribution coming from rotation is sufficiently large, the domain admissible for motion 
is empty. 

Finally, the smallest value of x — Xq is given by V(xq) — because the critical points must be always a zero of the 
potential function. 

From the physical point of view the critical point of this type represents the static Einstein universe. 
Let us briefly comment now on the case of vanishing rotation in the above context. Then the domain admissible 
for motion is for positive x. In the phase space all critical points are situated at the a>axis, at the point 



x 



Ofc,o 

3Oa. ' 



vo = o. 



Therefore, we have always static critical points if only £lk,o and Oa,o are of opposite signs. 

The character of the critical point is determined from eigenvalues of the linearization matrix at critical point or 
equivalently from the convectivity of V{x) 



d 2 V 
dx 2 



= -'SQ u ,qx 4 - fim.oa; 3 - Oa.o 



but from (|57|l we have at the critical point 



Then for any xq we have 



d 2 v 



dx 2 

d 2 v 



(x 0l Q) 



-£l m fiX 3 — 40a. o < or in general 



dx 2 



= det A = -A1A2 



(»o,Q) 



where Xi (i — 1, 2) are eigenvalues of the linearization matrix A of the system 



.4 





d 2 V 
dx 2 



Then eigenvalues are solutions of the characteristic equation A 2 — (tr A)X + det A = for the matrix A. 

Finally, we show that in the physical domain our system with positive curvature k > has at least one critical 
point which is represented on phase plane by the saddle point. 

Let us consider the phase portrait of the system with global rotation and dust matter on the phase plane (x, y). We 
assume the existence of the non- vanishing cosmological term because the value Oa,o 7^ is required by consistency of 
the model with global rotations with SNIa data. 

If we consider the standard FRW cosmology on the (x, y) plane then situation is presented in Fig. The flat 
model trajectory separates the region of the model with negative and positive curvature. The critical point is a 
saddle and it represents the static Einstein solution. Non-generic cases corresponding to separatices going in or out 
of the saddle point. The acceleration region is situated on the right from the saddle point. Therefore the recollapsing 
and then expanding models lie permanently in the acceleration domain and there is a transition from decelerating to 
accelerating epochs. The Eddington model is also in this region. The Lemaitre-Eddington models start accelerating in 
the middle of the quasi static phase. As we can see from Fig.^the global rotation contribution introduces qualitative 
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-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



x 

FIG. 1: The phase portrait of system for models with global rotation fi^o 7^ and dark matter X in the form of solid dark 
energy [vox = 1/3), i.e., x — y, y — — |fito P ,o — 30.^, ox~ 7 + Qa,ox where fito P ,o is density parameter for solid dark energy. 
The acceleration region is situated on the right from the saddle point. Therefore, the trajectories of recollapsing and then 
expanding models lie permanently in the acceleration domain. Algebraic curves, on which lie trajectories, are given in explicit 
form from of first integral ^- = §fifc ( o + ^flto P ,ox + ^Q^fiX -8 + ^Qa,ox 2 . The physical region is determined by the condition 
x > xq : V(xq) = where xo is the critical point. 



changes of behaviour of trajectories, i.e., there is no homeomorphism which preserves directions along trajectories 

and transforms the corresponding the trajectories. In our case there is no critical points on the phase plane at the 

1/3 

finite domain (Fig. It can be shown that in general if only a > (ipf) = a m in, then there is no critical point 
representing a static Einstein universe. 

To analyse system (I55|l at infinity it is useful to introduce the projective map (v,u): v = 1/x, u = y/x on the plane 
and 

v = —uv (59a) 
u = -7^m,o« 3 - ^,ov 4 + ^a,o - u 2 (59b) 

Then, circle at infinity is covered by line v = [x — 00). The dynamical system with the non-zero cosmological 
constant and non- vanishing global rotation in the finite region is presented in Fig. EI The region x < x m i n is forbidden 
because the global rotation contribution cannot dominate the matter contribution. In this case there is one type of 
trajectories which recollapses from the initial singularity and then expands after reaching the minimum size of scale 
factor. Hence the expanding de Sitter model is a global attractor and the contracting de Sitter one is a global repeller 
(Fig. |SJ. All point at infinity which lies on v = circle are hyperbolic, therefore the system is structurally stable. 
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FIG. 2: Diagram of the £Ia{x) density parameter given by relation: Qa{x) = — °' k -" J ' TJ '" 4 '" a ' for different values of (fifc.o, 

^m,o, fiu,o). Forbidden (non physical) region: y 2 /2 < is situated under the characteristic curve. The lines Q.a{cl) = const 
gives qualitative classification of possible path evolutions in the configuration space. The minimum is reached for Qk,o — at 



Note the existence of oscillating models without the singularity for £Ia,o < 0. 



VII. THE FLATNESS AND LAMBDA PROBLEMS 



If we choose the standard FRW cosmology without the global rotation then the conservation equation for perfect 



fluid ©, satisfying the equation of state pi = (7$ — l)pi, gives 



Pi oc a 



-3 7 » 



and the p/3 term, in the Friedmann first 



integral (called also the Friedmann equation). Let us consider first these problems by setting A = 0, 



= -2V(a) = f 



2 to^- w 

3 (3 7 - 4) 



k 



(60) 



where ao = 0, i.e., shear vanishes. 

The matter term dominates the curvature term ka~ 2 at large as long as the matter stress obeys p + 3p < 0, 
P + P > 0, that is ifO < 7 < 2/3. This is what we describe as the flatness problem. Since the scale factor then 
evolves as a{t) oc t 2 ^ 31 if 7 > or if 7 = 1 as exp(H t). We can observe that it grows faster than the proper size of 
the particle horizon scale (oc i) so long as 7 < 2/3. Therefore, a sufficiently long period of evolution during which 
expansion is dominated by matter with negative pressure with p + 3p < can solve the flatness and horizon problems. 
Such a period of accelerated expansion is called inflation 43] . However, as Barrow noted there is a flatness problem 
if the cosmological constant is added to the right hand side of the Friedmann equation. Then, to explain why it does 
not dominate matter term p/3 (at large a(t)) it is assumed that there exists the period in the early evolution of the 
universe during which weak energy condition for the matter is broken, i.e, p + p < 0. This is what Barrow called the 
cosmological constant problem. For the standard FRW model the existence of corresponding 7 < from equations 
of state mean that a < which is the case of phantoms. Therefore, models with global rotation provide a solution of 

the flatness problem. It is formally equivalent to presence in the model some fictitious additional noninteracting fluid 

67-10 



which obeys the equation of state p^ — 7 3 67 p, 



(a ) 



Pu a — 



37^4^0 ■ 



Then both weak and strong 




2 4 6 8 10 

x 



FIG. 3: The phase portrait of the FRW model with global rotation, the cosmological constant and dust given by system 
All trajectories start from state x = oo,x — — oo and after bouncing expand to infinity. Cosmologies with rotation are 
represented only by physical part of the phase curve starting from point x = x m i n ,y = 0. There is no critical point in the 
physical domain. 



energy conditions for pure "rotational fluid" means 

IO-67 

Pu + Pu = w u p u = p u > 

A, + 3^ = 2/^(4- 3 7 ) >0. 
In the special case of dust matter we have inequality if the above conditions are violated 

4 -4 
Plo + Puj = T^Puj < 0, Pu=Pu X 

Plo + ZPuj = ^Plo < 0, p Wo = -2u%, 

i.e., both energy conditions are violated. 

Finally, the strong and weak energy conditions for density energy of dust and rotational fluids requires respectively 

(Pm + pJ) + 3p w = Pm,oa -3 (l — ujIcl^ 1 ] > 0, 

V Pm,0 / 
(pm + Pu) + Pu = Pm.oa~ 3 (l - T~~ UJ O a l ^ °- 

Thus, during the evolution, first, the strong energy condition is violated for a < -^-^ujq = n 1 - and the weak 
one is violated for a < 3p s - . We can also find the interval 2p < a < 3p s ^ luq in which the strong energy condition 
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FIG. 4: The phase portrait of system 155H with the cosmological constant but without global rotation f2 w ,o = 0. The acceleration 
region is situated on the left from the saddle point. Therefore, the recollapsing and then expanding models lie permanently 
in the acceleration domain. The Eddington model also belongs to this region. The Lemaitre-Eddington type models with 
characteristic loitering phase start accelerating in the middle of quasi-static phase. 

is violated, whereas weak energy condition is still satisfied. It is possible because p w is negative. These are sufficient 
conditions to solve the flatness problem. 

Because of the existence of the period in which p + p < (phantom), the cosmological constant problem can be 
solved and a < is an unnecessary requirement. Unfortunately, the future horizon problem cannot be solved and the 
corresponding condition 7 < 2/3 is identical to that of the solution of the horizon (or flatness) problem in the case 
with vanishing rotation. 

Let us note that if we consider the case of the non-zero A term, that is contributed by a pressure p\ = — pa then 
A term on the right-hand side of (|60[l falls off faster than the curvature and matter density terms so long as 7 < 2/3 
and 7 < respectively. The contribution term from the global rotation falls off faster then curvature A and matter 
density term so long as 7 < 4/3 (radiation), 7 < 5/3, 7 < 10/9, respectively. 

Therefore, the cosmological constant term A, curvature term is proportional to a~ 2 and term connected with global 
rotation falls off faster than 8wGp/3 term in the FRW equation at large a(t) if respectively 7 < and 7 < 2/3 and 
7 < 10/9, i.e., 7 < 0. 



VIII. HOMOLOGOUS UNIVERSE AND ITS ANGULAR MOMENTUM 



In this section we demonstrate how Wesson's argument can be generalised and how it works in our case |48l |49j . 
Wesson argued that because the dynamics of a gravitationally bound rotating system admits self similarity group 
of symmetry the relation between angular momentum J of a astronomical system and its mass J ~ A/ 2 can be 
constructed due to constancy of the invariants of this kind of symmetry. In this section we demonstrate how Wesson's 
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FIG. 5: The phase portrait of the model 15911 with global rotation fi Wj o 7^ 0, the cosmological constant and dust in the projective 
coordinates (v,u): v = 1/x, u — y/x. By introducing such coordinates we can observe their behaviour at infinity. There are 
two types of trajectories which expand to the maximum scale factor and then recollapse and vice versa. In the first type of 
evolution the closed models start with the singularity a(0) = and a = 0, then reach the maximum radius and recollapse while 
the other class of models recollapse from a(0) = oo and a — — oo to the minimum value of a, and then expand to infinity. In 
the generic case trajectories start from the anti-de Sitter stage and land at the de Sitter state. 



argument can be generalised and how it works in our case. We generalise Wesson's argumentation to the case of 
universe dynamics by consideration of Lie (continuous) symmetries of the FKW models with rotation. The Lie group 
theory of symmetries analysis gives us information about the group of symmetry transformations admissible by the 
differential equation structure. They preserve the structure of basic dynamical equations. In the set of solutions, 
the action of an admissible group includes a certain algebraic structure which can be used to find a family of new 
solutions from the known ones [44L \43L I46L l47j . 

As is well known [50| , the new solution for the system equation, describing the static stars (in Newtonian as well as 
in General Relativity), can be obtained from the known ones through the homologous transformation |5lll52ll53| . We 
can find the homologous transformation of the symmetry type (i.e., preserving structure of differential equations) for 
the considered dynamics of the FRW model with rotation. It would be useful to rewrite basic dynamical equations 
l{9"jl and J3J to the form in which instead of (H, p, u>) variables we have (a, p, u). Then we obtain 



dt V 3 3 v ' 



± = -2<yp\ P —-k- -^a^a- 1 (61b) 
dt V 3 3 y ' 



^ = -o;(5 - 37)^ - A - \^a~\ (61c) 
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Now let us consider the differential equations system 

— = f(x,u ,...,u m ), i=l,...,m, 

and here the space of independent variable is denoted as x, dependent one u a and its first derivatives, say u a (a — 
1, . . . , to). The action of the Lie group G of point transformation in the space {x, u) is described by an infinitesimal 
operator X — a generator of the symmetry 

a TO o 

X = C(x,u\...,u^—+Y,v\x,u\...,u m ) — . 

i—l 

The point transformation generated by X is called homologous if Q — ax and rf — g l u l , where a,g l (i = 1, . . . , m) 
are constants. Of course, if X = X s (u)d s , then the finite transformation of symmetry u — » u are given as a solution 
of equation = A s (u) with the initial conditions u s (t = 0) = u s (s = 1, . . . ,m). The action of the Lie group 
transformation can be extended from the (x, u) space to the (x, u, u') space, i.e for the first derivatives. On the other 
hand, s-th order (ordinary or partial) differential equation 

F(x, u{x), u l {x), . . . , u s (x)) = 

defines a certain manifold M. in the space (x^UjU 1 , . . . ,u s ). We say that F is invariant with respect to the action 
group G, provided that manifold M. is a fixed point with respect to the s-th extension of G, i.e., G S (A4) = A4. In 
the terms of the infinitesimal operator it means that 

X S F\ F=0 = 

where X s is the extented operator on the first derivatives 

X^X + Q{ x^u>) d — f < = g, 

Cf = A(rf) - upX^)., D l= d —+ uf—. 

The prolonged operators, which are generators of symmetry of equations in the space (x, u), form the structure of a 
Lie algebra of the fundamental group. 

The Lie method adopted to our case gives us the symmetry operator in the form 

d a d a d a d , 
X = a %-2 a d^ + 2 UJ d^-2 t d-t' (62) 

where a = const, i.e., the most general form of symmetries constitutes a homologous transformation. For our further 
analysis it is useful to introduce the notion of invariants of symmetry operator X. For the infinitesimal operator X 
there is to independent invariants which are solutions of the following system 



dx du du T ' 



(63) 



Q{x, u 1 , . . . , u m ) rfix, u 1 , . . . , u m ) r) m (x, u 1 , . . . , u m ) 

In our case we can recover two invariants in the form 

J\ = pa 2 , J2 = toa. (64) 

They exactly correspond to Wesson's invariants expressed in dimensionlcss form 771 = Gpa 2 /c 2 , r/2 = uja/c. The 
significance of the obtained invariant homologous transformation of l|45|) is the following. 

1) They can be used as good new variables to integrate the system by lowering its dimension. Then we obtain 

^ = (2-37)Ji, ^ = (3 7 -4)J 2 , r,:Hdt = d V . (65) 
dr\ dr\ 

2) The invariant can be used to obtain a new solution from the known ones using the homologous theorem. If a(t) 
is a solution of (|45|1 . then because a(i)i _1 is invariant, a(e~ aT / 2 t)(e~ aT / 2 t)~ 1 = const is also a solution, where the 
finite transformations are given in the form 

p ^p = pe <xT : a ^a = ae- aT ^ 2 , uj -> u = ue aT ^ 2 , t^i=e~ aT / 2 (66) 
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where r is a group parameter and a = const. Let us note that the finite transformation (|66|l preserves the angular 
momentum form 

J = pa 5 oj — pa 5 cD, (67) 

i.e., this relation is invariant with respect to the homologous (or self similarity) transformation group given by i|66[l . 
Physically the existence of this type of symmetry always means that a basic dynamical equation can be formulated 
in a dimensionless form. 

It is interesting that (|(jT|) can be expressed in terms of invariants by 

J=^-M 2 , (68) 
Jl 

and astrophysical relation J ~ M 2 which is closely obeyed over mass range 10 18 <j to 10 48 <7 can be extended to 
the Universe, together with the Wesson's philosophy that self similarity transformation plays an important role in 
astrophysics. This notion was introduced in astrophysics by Rudzki |54j. We note that it can play an important role 
in cosmology (see also 0])- 

Moreover the ratio Jij J\ in the equation (|68|l is constant in the distinguished case of dust matter (7 = 1). 
Therefore, there is no need for additional arguments which usually are taken from the astronomical observation about 
the constancy J2/J1 as in the Wesson argumentation because after integration of (I65|) we obtain immediately 

J2 



Ji 

and only for 7 = 1 (dust) we obtain 

J2 



e 



6(1-7)*? 



(69) 



T = const. (70) 
Ji 



IX. CONCLUSION 



In this paper we studied the dynamics of the FRW universe with global rotation. Our approach is simplest in that 
we formulate the dynamical problem in two-dimensional phase space. Moreover, we find the Hamiltonian formulation 
of the dynamics. Such visualisation has a great advantage because it allows us to analyse the acceleration problem in a 
clear way. Due to the existence of the Hamiltonian constraint, it is possible to make the classification of the qualitative 
evolution paths by analyzing the characteristic curve which represents the boundary equation in configuration space. 
Dynamical system methods gives us all evolutional paths for all possible initial conditions. On the other hand, the 
representation of dynamics as a one dimensional Hamiltonian flow allows us to make the classification of the possible 
evolution paths in the configuration space. When we consider the dynamics of the FRW models with global rotation, 
a two-dimensional dynamical system in general, then there is a simple test of structural stability of the system (a 
dynamical system S is said to be structurally stable if dynamical systems in the space of all dynamical system which 
are close to S are topologically equivalent (physically realistic models on the plane should be structurally stable). 
Namely, if the right-hand sides of dynamical systems are in polynomial form, the global phase portraits are structurally 
stable on S 2 {R 2 with adjoint a Poincare sphere) if, and only if, the number of critical points and limit cycles is finite, 
each point is hyperbolic and there are no trajectories connecting saddle points. Our conclusion is that the FRW 
models with global rotation are structurally stable, i.e., following Peixoto's theorem dynamical systems with rotation 
form open and dense subsets in the space of all dynamical system on the plane. From the point of view of science 
modelling they are better and more adequate description of the observed universe than the FRW model without the 
cosmological term. 

Our formalism gives natural base to express dynamical equations in the form of the FRW model with some addi- 
tional fictitious noninteracting multifluid which mimics the effects of global rotation. It satisfies in the generic case 
the equation of state for radiation with some unusual negative energy density. We find basic dynamical equations 
describing the dynamics in the form of two-dimensional Hamiltonian dynamical systems in which coefficients are just 
dimensionless observational density parameters Q^o- 

We also showed that global rotation produces a dark energy component but the cosmological term is still required 
to explain SNIa data. We also demonstrates how classical cosmological problems can be solved due to the existence 
of global rotation. We also derive the observationally suggested relation between the angular momentum and mass 
of the object J oc M 2 from the property of self similarity of dynamics of the FRW cosmology with rotation. 

The presented formalism forms the base for our further analysis of the magnitude redshift relation and finding the 
best fitting for f2 Wj o from recently available SNIa data and X-ray gas measurements of mass fraction in clusters of 
galaxies or measurements of angular size of radiogalaxies. 
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